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ABSTRACT 

Cosmological transverse momentum fields, whose directions are perpendicular to Fourier wave vec¬ 
tors, induce temperature anisotropies in the cosmic microwave background via the kinetic Sunyaev- 
Zel’dovich (kSZ) effect. The transverse momentum power spectrum contains the four-point function 
of density and velocity fields, {6Svv). In the post-reionization epoch, nonlinear effects dominate in the 
power spectrum. We use perturbation theory and cosmological A^-body simulations to calculate this 
nonlinearity. We derive the next-to-leading order expression for the power spectrum with a particular 
emphasis on the connected term that has been ignored in the literature. While the contribution from 
the connected term on small scales {k > 0.1 h Mpc“^) is subdominant relative to the unconnected 
term, we find that its contribution to the kSZ power spectrum at £ = 3000 at z < 6 can be as 
large as ten percent of the unconnected term, which would reduce the allowed contribution from the 
reionization epoch (z > 6) by twenty percent. The power spectrum of transverse momentum on large 
scales is expected to scale as as a consequence of momentum conservation. We show that both the 
leading and the next-to-leading order terms satisfy this scaling. In particular, we find that both of 
the unconnected and connected terms are necessary to reproduce k^. 


1. INTRODUCTION 


The cosmic transverse momentum field is observable in 
temperature fluctuations of the cosmic microwave back¬ 
ground (CMB). The line-of-sight (LOS) component of 
the momentum of electrons in the ionized intergalactic 
medium (IGM) induces temperature fluctuations via the 
Doppler effect, and this i s known as the kinetic Sunyaev - 
Zel’dovich (kSZ) effect (|Zerdovich and SunvaevI 119691) . 
The contribution from the longitudinal momentum field, 
whose direction is parallel to the wave vector, suffers from 
cancellation of the contributions of troughs and crests in¬ 
tegrated along the LOS; thus, the cont ribution from t he 
transverse momentum field dom inates (iVishn iadllQS Ttl. 

Th e kSZ effect is given by (iSunvaev and Zerdovichl 

Fl^ 


T J c ’ 


( 1 ) 


where 7 is the LOS unit vector, v the peculiar velocity 
field, and r the optical depth to Thomson scattering in¬ 
tegrated through the IGM from z = 0 to the scatterer. 
The differential form of r is dr = c necrx (^) dz, where 
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rie is the number density of free electrons in the IGM, 
and ax is the Thomson scattering cross section. 
Equation © can be rewritten in the following form: 



arUefi f ds 

- / — e q -7, 
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( 2 ) 


where q = xv(l -|- <5) is the momentum of ionized gas, <5 
is the fractional mass density perturbation of the gas, 
X = ne/[n-ii + 27T,He) the ionization fraction, he,o = 
fi-H.o + 2hHe,o the mean number density of electrons at 
the (fully-ionized) present epoch, and s the distance pho¬ 
tons travelled from a source to the observer in comoving 
units. 

Longitudinal momentum fields cancel out in the line- 
of-sight integral of Equation ([2]) (See Appendix lAl for a 
quantitative argument; also see lVishniaalI987tl . Then, 
the angular power spectrum of Equation Q, at 
large multipoles is dominated by the power spectrum of 
the transverse momentum field, {k). and is given by 
dVishuiaHhflSTh 
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The kSZ power spectrum has been constrained by ob¬ 
servations. The South Pole Telescope (SPT) and Ata¬ 
cama Cosmology Telescope (ACT) measure CMB tem¬ 
perature anisotropy at scales beyond the Silk damping 
scale {£ > 3000) where it is possible to distinguish the 
kSZ signal from the primary signal. The latest con¬ 
straint is = 2.9 ± 1.3 uN (6 9% CL), where 

Dg = l{i + l)C'^/27r (jCeorge et alll20lB l. 

The measured kSZ power spectrum is the sum of con¬ 
tributions from the epoch of reionizat ion (EoR) and the 
post-reionization epoch. The model of iShaw et al\ (l2012f ) 
suggests that the latter gives ~ 2.0 as¬ 

suming that the universe became fully ionized at z = 6. 
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Thus, the central value of the current measurement im¬ 
plies that the post-reionization kSZ signal is a factor of 
two greater than the EoR signal. This highlights the im¬ 
portance of understanding the post-reionization signal; if 
we mis-calculate the post-reionization signal by ten per¬ 
cent, the EoR signal would be mis-estimated by twenty 
percent. This motivates our revisiting assumptions used 
by the previous calculations of the post-reionization kSZ 
power spectrum. 

Inhomogeneity in the ionization fraction during the 
EoR gives a large boost in the kS Z power spectrum com¬ 
pared to homogeneous ionization (jPark et a/.ll201,1[ l. The 
physics that determines this inhomogeneity is complex, 
and many papers have been written on this subject, rang¬ 
ing from early analytica l calculations l)Gruzinov and Hul 
llQQStl^ntos et alll2003[l and semi-numerical calculations 
using an analytica l ansatz applied to A^-body simulations 
(iZahn et oi.ll200,^ iMcQuinn et oi.ll2005t iMesinger et al\ 

l2012HZahn et al.ll2011l : lBattaglia et al.ll2013ll . to fully nu- 

merical simulations that couple the co smological struc¬ 
ture formation w ith radiative transfer (|Iliev et al.ll2007t 
iPark et "^1201,il l. 

Modeling the post-reionization signal is simpler be¬ 
cause the IGM is fully ionized and y does not fluctuate 
to a good approximation. We thus need to model the 
density and velocity fluctuations of gas. When y = 1, 
the momentum held is given by q = v(l -|-(5). Since 
the velocity held is purely longitudinal in the linear 
regime, the transverse momentum held, q_L, is given 
by q± = (v(5)_l at leading order. The power spec¬ 
trum is then given by the four-point function of two 
densities and velocities. Schematically, it is given by 
( 99 ) = {vv){SS) + 2{v5)^ + {vSvS)c- The last term is 
called the connected four-point term, while the others 
are the unconnected ones. 

In the previous work, the connected term has been 
ignored. For example, the earlier analytical studies 
ignore nonlinearity in density or velocity, and calculate 
only the unconnected terms u s ing li n ear perturbation 
theory (lOstrike r and Vishniad [198111 IVishniad Il987t 

Dodeison and Jubasl 119951 I.Taffe an^ <amionkowskil 

1998f) . An analytical model of IHiJ ( 200011 is still 
based upon linear theory for velocity, and ignores the 
connected term, but replaces the linear density power 
spectrum wit h a model for the nonlinear d ensity power 
spectim m by iPeacock and Dod^ (|1996l l. iMa and FrvI 
(I2?i0^ use a similar approach with a halo model for 
the nonlinear density power spectrum and argue that 
the connected term is negligible at large k. While we 
broadly agree with this conclusion, our aim is to quantify 
the contribution of the connected term at large k, and 
also clarify the role of the connected term in obtaining 
the correct small-A: limit of the transverse momentum 
power spectrum in perturbation theory. 

Some of the previous “numerical” calcul ations of the 
post-reionization kSZ power spectrum (I Zhang et al\ 
12004 iShaw et o/.ll2012ll still rely on the above analytical 
model that ignores the connected term, but takes the 
ingredient of the model, i.e., nonlinear gas density power 
spectrum, from hydrodynamical simulations. Therefore, 
quantifying better the contribution of the connected 
term affects t he results from the previ ous numerical 
work as well. iSoringel et al\ l)2001h and Ida Silva et al\ 
(|2001ll computed the kSZ power spectrum directly from 


their simulation and did not rely on the model. 

Throughout this paper, we shall assume that gas traces 
dark matter. This is not a great approximation; shocks 
in the IGM generated b y structure formation he at gas to 
high temperatures le.g.. lGen and Ostrikerlll99^ . and the 
resulting gas pressure makes gas less clustered than dark 
matter particles. This effect on the kSZ power spe ctrum 
is modest at ^ = 3000 (|Shaw et aLll2012l : lHull200dll . Star 
formation converts gas into stars, further reducing the 
kSZ effect. Shaw et al. (2012) find that can be 

suppressed by as much as 33% of the simulation without 
gas cooling and star formation. Our goal in this paper is 
to quantify the error we make by ignoring the connected 
four-point term in the transverse momentum power spec¬ 
trum. While our dark-matter-only results cannot be ex¬ 
trapolated to gas, we expect that a similar conclusion 
would still apply to gas, at least qualitatively. 

The remainder of this paper is organized as follows. In 
Section[21 we briefly introduce the A^-body simulation we 
use in this paper. In Section [3l we review the derivation 
of the current n onlinear transverse momentu m power 
spectrum model (jHd 1200(11: IMa and Frvll2002h . In Sec¬ 
tion |4j we first confirm that the current model accurately 
approximates the unconnected term. We then show that 
our simulation data of the transverse momentum power 
spectrum exceed the model. We argue that this excess 
is due to the connected term, by showing that the per¬ 
turbation theory calculation of the next-to-leading order 
terms (including the connected term) explains the excess 
at quasi-linear scales successfully. We also show that the 
connected term is essential in obtaining the correct low- 
k limit of the transverse momentum power spectrum in 
perturbation theory. In Section [SJ we quantify the im¬ 
pact of the connected term on the post-reionization kSZ 
power spectrum. We conclude in Section HI 

2. SIMULATION 

We shall use a cosmological 7V-body simulation of col- 
lisionless particles using the “CubeP^M” 7V-body code 
(jHarnois-Deraps et a/.l[2f)T^ . The simulation is run with 
3456^ particles in a comoving box of 1 h~^ Gpc on a 
side and is started at 2 = 150 using the initial condition 
generated using the Zel’dovich approximation and initial 
density power spectrum fro m the publicly available code 
GAME (iLewis et al.l 1200011 . Thi s simulation was previ¬ 
ously presented in iWatson et al\ (|2013f) . 

The resolution of this simulation allows us to sample 
on average ^ 40 particles per (Mpc//i)^. This resolu¬ 
tion allows us to avoid sampling artifacts in the velocity 
power spectrum u p to /c ~ 1 /i Mpc“^ (see Figure 3 of 
iZhang et al.ll2015h . We then adaptively smooth particles 
to a grid of (1008)^ cells 0. Therefore, our simulation 
covers a dynamic range of 0.006 ^ k < 1 h Mpc“^ in 
wavenumber. The background cosmology is based on 
the WMAP 5-year data combined with constraints from 
baryonic acoustic oscillations, from observations of galax¬ 
ies and large-scale structures, and from high-redshift 

® This adaptive smoothing is by an approach similar to that 
used in Smoothed Particle Hydrodynamics. In this case, spher¬ 
ical smoothing kernels are assigned to each particle, with their 
smoothing lengths adjusted so as to enclose the locations of the 
32 nearest-neighbor particles. The mass per particle assigned to a 
given grid cell then corresponds to the fraction of its finite kernel 
volume which overlaps the cell volume. 
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Type la supernovae {flM = 0.279, r^A = 0.721, fe = 
0.701 , flh = 0.0462,(78 = 0.8, rig = 0.96: iKomatsu et~^ 
[200^ . 


3. TRANSVERSE MOMENTUM POWER SPECTRUM 
3.1. Basics 

In the post-reionization era, helium atoms are singly 
ionized until He II reionization occurs. We assume that 
hydrogen reionization finished at z = 6 and He H reion¬ 
ization occurred instantaneously at z = 3; thus, x = 0.93 
for 3 < z < 6 and y = I for z < 3. For the rest of the pa¬ 
per, we shall drop x foi' notational simplicity and write 
q = v(l -|- 5). Therefore, the momentum power spec¬ 
trum below should be rescaled by x^ to yield the correct 
value needed for computing the kSZ power spectrum in 
Section [S] 

We start by Fourier transforming the momentum field, 
q = v(I -I- (5), as 

q(k) = v(k) + J j^S{k - k')v(k'), (4) 

where q(k) = J e*'^'^q(x), etc. Then, its power spec¬ 
trum is defined by 


(2^)3p*^(ki)d,3(ki +k2) = (gXki)g'(k2)) 

^(d(ki-k'y(k')u^(k2) 


d^k' f d^k" 


(277)3 


(277)3 


d(ki - k')ii*(k')d(k 2 - k!')v^{k") 


(5) 


where the indices, i and j, denote the I’th and j’th com¬ 
ponents of the vector, respectively. The power spectrum 
of the transverse mode, qj_ = q — k[q • k], is given by 
Pq^ = [1 —(fc*)^]. In linear and quasi-linear regimes 

(fc < I h/Mpc), the velocity field is longitudinal (i.e. 
v(k) = ku) to a good approximation, which implies that 
{vv) and {5vv) in Equation ([5|) vanish. Figure [T] shows 
that the contribution to Pq^ from {5vv) (diamonds) is a 
few orders of magnitude smaller than {SSvv) (triangles). 
The contribution from (vv) (not shown) is even smaller 
than {Svv) by another two orders of magnitude. 

Keeping only the last term in Equation ([S|) , we obtain 


z = 0 



Fig. 1.— Contributions to k from (diamonds) and 

(triangles) in Equation (JS]! measured from our simulation 
at 2 = 0. The red triangles show negative values. The dashed line 
shows the lowest order calculation for the {(55r;r;) t erm , while the 
dotted line shows its low-A: limit given in Equation {nj). 


= (5(ki - k')u*(k')) (^(k2 - k'Ou^k")) 

+ (5(ki - k')5(k2 - k")) (u*(k')u^ (k")) 

+ (j(ki - k')ii^(k")) (5(k2 - k")u*(k')) . (8) 


The firs t term in the abov e vanishes and the other terms 
lead to (|Ma and Ervll200^ 


P,Ak,z)= J [;^P^^(|k-k'|)P,,(fc') 


Ik-k 


-^Pse{\k-k'\)Pse{k') 


, (9) 


where 0 = V • v. 

The linear velocity is related to the linear density 
hy 9 = df6, where / = din 6/din a and a(t) is the 
Robertson-Walker scale factor. This gi ves the lowest or¬ 
der (LO) expression as (|Vishniaclll987ll 


/ d^k' f d^k" ~ 

P55™(k-k',-k-k",k',k"), (6) 

where Pssw is defined by 

(27r)3p55„„(ki,k2,k3,k4)5£,(ki + k2 + ks -b k4) 

= (l(ki)5(k2)i}*(k3)h^(k4)). (7) 


3.2. Linear Regime 

In the linear regime, we consider only the first order 
terms of ^’s and u’s in {66vv). Then, Gaussianity of linear 
V and S fields yields 

(5(ki - k')u*(k')^(k2 - k")i)^(k")) 



fc(fc-2fcV)(l-/i'^) 
fc'2(fc2-bfc'^-2fcfcV)’ ^ 


where is the linear matter density power spectrum. 
A similar derivation for the power spectrum of the sum 
of longit udinal and t ransverse momentum fields is pre¬ 
sented in iParkI (|200Clll . 

Taking /c —^ 0 limit of Equation we find 


P^°ik,z) 


15 




dk' [P^i^{k',z)f 

(277)2 -11 


The k dependence of Pj(® is thus given simply by 
which is independent of cosmology or the initial power 
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spectrum. Similarly, the second order term (2nd term in 
the r.h.s. of Eq. lAlOl) of the longitudinal mode also goes 
as in the low-Zc limit: 




dk' [P^Tik',z)? 

(27r)2 fc/2 • 


We shall discuss the physical origin of this dependence 
in Section 

We use CAMB to compute The dashed line in 

Figure [1] shows at z = 0, while the dotted line shows 

the low-fc limit lEg. fTTI) . 


3.3. Nonlinear Regime 

As Figure [1] shows, nonlinear contributions to 
Pq^ become dominant over the LO contribution at 
k > 0. 2 h/Mpc at z = 0, in agreem ent with the previous 
work (|Hull2000l : iMa and Frvll2002ll . The current popular 
model of the post-reionization kSZ signal uses an 
approximate expression for the t ransverse momentum 
power spectrum due to Ih71 (HooO), which replaces one 
of in Equation (TTOl) with the nonlinear power 

spectrum, P^J: 

plik,z)=d^f j - k'l,z)Pir(fc', 

fc(fc-2fcV0(l-/i'^) 

P'^ik^ + kP -2kk’^i')' ^ ’ 

This model modifies the second term in the square 
bracket of Equation ((O]) in the following way: 

Pse{\k - k'\)Pse{k’) = Pss{\k - k'\)Pss{k'). ( 14 ) 

This holds in the linear regime, and the entire term 
vanishes in the large k regime due to the pre-factor, 
|k — k'|“^. In addition, the model approximates the 
velocity power spectrum by linear theory, i.e., Pgg = 
a? pPgP ■ In this way, the model avoids having to model 
the velocity power spectrum and the density-velocity 
cross power spectrum, which is relatively poorly under¬ 
stood. We shall refer to this model as the Standard model 
(hence the superscript “S” in Eq. [T5]) in this paper. 

However, Equation o, which the Standard model 
aims to model, is not the full expression for Pq^^ be¬ 
cause it neglects the contribution from the connected 
term, {SvSv)c- In the nonlinear regime, nonlinear growth 
makes both S and v non-Gaussian, and thus there is no 
reason to think that the connected four-point term is neg¬ 
ligible. We shall quantify the importance of this term in 
detail in Section 021 

4. REVISITING THE STANDARD MODEL FOR THE 
EFFECT OF NONLINEARITY ON THE TRANSVERSE 
MOMENTUM POWER SPECTRUM 

4.1. Uneonneeted Term 

In this section, we revisit accuracy of Equation (IT3l) . 
As noted above, this model is an approximation for the 
unconnected term (Eq. [9l). Therefore, it, by design, does 
not take into account the connected term. 

We test whether Equation (021) successfully approxi¬ 
mates Equation by evaluating it using the density 
power spectrum from the simulation, i.e., P^^ = Pgg^, 


and compare it with Equation ([S]) using Pss, Psg, and 
Pgg from the simulation. 

In Figure [U we show the ratios of Equation m 
(squares) and Equation ([9]) (crosses) to Pqj^ measured 
directly from the simulation. At 0.2 ^ k < 0.5 h Mpc“^ 
at z = 0 (top left panel), the Standard model repro¬ 
duces the unconnected term with high accuracy (~ 2% 
level). However, the Standard model overestimates the 
unconnected term at larger k, reaching ^ 5% level at 
fc ~ I h Mpc“^. We attribute this error to the linear ve¬ 
locity assumption overestimating magnitudes of velocity 
modes. I n the high-fc limit, Equation ([9|) converges to 
mull2000ll 

PqAk) = \Ps&(k)vlms(k), (15) 

where vl^^{k) = Jk><k-0^Pvv{k') is the velocity 

dispersion. Nonlinear correction makes the velocity 
power spectrum in the relevant k range smaller than 
the linear velocity power spectrum (see, e.g., Figure 
I of iPueblas and Scoccimarrol 12009^ : thus, linear the¬ 
ory overestimates This effect becomes smaller 

at higher redshifts, as shown in the other panels of 

Figure [21 

We find that both Equation ([2]) and uni) underesti¬ 
mate Pqj^ significantly compared to Pqr measured di¬ 
rectly from the simulation. The underestimation de¬ 
creases monotonically with redshift: 12 — 25% over k = 
0.2 - 1.5 h Mpc"^ at z = 0; 7 - 23% at z = I; 5 - 18% 
at z = 2; and a few to 5% at z = 5.5. Thus, this is 
likely related to the development of nonlinear structure 
formation. 


4.2. Connected Term 

To show that the connected term is the likely expla¬ 
nation for the difference between the simulation data 
and the models that ignore the connected term, we cal¬ 
culate the connected term using perturbation theory. 
Since the connected term vanishes in linear theory, we 
must go to the next-to-leading order perturbation the¬ 
ory, such as the standard “o ne-loop” perturbation the¬ 
ory ([Bernardeau et a?] 1200211 . This theory allows us to 
extend validity of analytic solutions for Fourier modes 
of density and velocity fields down to weakly nonlinear 
scales, e.g., in /c < 0.3 h Mpc“^ at z ^ I and in wider 
wave numbers at higher redshifts (iJeong and Komatsul 
1200611 . We also confirm that the density power spectrum 
from our simulation data supports their findings. We de- 
rive the explicit expressions for t he connected (Eq. IBI6|) 
and unconnected terms (Eq. IC2[) in Appendix 0 andICl 
and show them in the dotted and dashed lines in Fig¬ 
ure 01 respectively. 

We find that each of the unconnected and connected 
terms does not vanish in the low-A: limit on its own. In¬ 
stead, they converge to a constant value with opposite 
signs. This const ant ca n be obtained by taking /c —>■ 0 
limit of Equation (IBI6I1 : 

„ , 16 f d^k' f d^k" [1-{k-k'){k-k")] 

qr,c\k^o-^ J J 

X P(k')P(k' + k")P(k") .(16) 
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Ratio to at z = 0 



Ratio to Pq.perpQt Z = 1 



k [h Mpc ’] 


Fig. 2.— Ratios of various models for Pqj^ to Pqj^ directly measured from the simulation at 2 = 0 (upper left), 2 = 1 (upper right), 
2 = 2 (lower left) and 2 = 5.5 (lower right). The models include the Standard expression (Ea. USD with P^^ from the simulation (squares), 
the unconnected term (Eq. EJ with Pss, Pse and Peg from the simulation (crosses), and the connected term calculated from perturbation 
theory added to the crosses (triangles). The vertical dotted line den otes the wavenumber beyond which the perturbation theory fails to 
model nonlinearity precisely according to IJeong and Komatsul 112009 ) and our data. At z = 5.5, that wavenumber is out of the plotting 
range and not shown in the figure. 


The sum of the two terms, however, yields a physical 
result that remains positive and decays toward lower k. 
The cancellation of two large numbers introduces an un¬ 
certainty in our numerical integration using the Monte 
Carlo methotfl. Within the uncertainty, the result is con¬ 
sistent with scaling, which is the same as the scaling 
of the LO expression given by Equation m- This calls 
for a physical explanation; namely, what is the physical 
origin of the scaling in the low-fc limit of the trans¬ 
verse momentum power spectrum, which is independent 
of cosmology o r the initial power spectru m? 

As shown by iMercolli and PaiCT ()2014[) , this is a con¬ 
sequence of momentum conservation. In short, the ar¬ 
gument goes as follows. Suppose that we have a uni¬ 
form distribution of matter particles with no initial mo¬ 
mentum or density fluctuation, and displace these par¬ 
ticles with velocities (where i is a particle ID) in a 

® In principle we can derive the low-fc limit of the next-to-leading 
order expression analytically. Here, we have chosen to perform 
numerical inte grati on be caus e of the complexity of the results given 
in Equations 1 IBI 6 I 1 and 


momentum-conserving manner. In this way, we have re¬ 
moved the effect of the initial condition, and can focus 
only on the effect of the subsequent evolution of particle’s 
motion. Fourier transform of momentum of the displaced 
particles is given by q(k) = exp(— zk • and 

the low-fc limit is q(k) = WiVi(I — zk • x,; -b ...). The 
first term vanishes by momentum conservation, giving 
q(k) = 0{k)\ hence the power spectrum of momentum 
fields is proportional to fc^. This argument applies to 
both the longitudinal and transverse momentum fields. 
As we have found, the connected term is necessary for 
obtaining the correct low-fc limit at the next-to-leading 
order in perturbation theory. 

5. IMPLICATION FOR THE POST-REIONIZATION KSZ 
ANGULAR POWER SPECTRUM 

In the high-fc regime, the connected term is of order 
ten percent of the unconnected term, which brings the 
predicted into better agreement with the simula¬ 
tion, as shown by the triangles in Figure [2j Thus, the 
underestimation is now much reduced to 2 — 18% over 
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0.0001 0.0010 0.0100 0.1000 1.0000 10.0000 
k [h/Mpc] 


Fig. 3. — Pqj^ at z = 0 from the next-to-leading order un¬ 
connected term (black dashed), connected term (black dotted), 
the sum of the two next-to-leading order terms (blue dashed), 
and the total sum of the leading and next-to-leading order terms 
(blue solid). The red color shows negative values. The er¬ 
ror bars show the uncertainty of the Monte-Carlo integration for 
k < 0.01 h Mpc“^. The solid line shows the LO term (Eq. llOH . 
while the blue dotted line shows its \ow-k limit fEo. Ill|l . 


k = 0.2 — 1.5 h Mpc“^ at z = 0; 2 — 18% at z = 1; 
2 — 14% at z = 2; and a few percent at z = 5.5. For each 
redshift, we mark in the figure roughly the wavenumber 
that the next-to-leading order (“one-loop”) calculation 
of perturbation theory beco mes inaccurate for the den - 
sity fluctuations according to lJeong and Komatsul 1)20061) 
and our data. The remaining differences beyond that 
wavenumber are likely due to inaccuracy of perturbation 
theory calculation, whereas the agreement at 2% level 
at small k supports our conclusion that the connected 
term is necessary for accurate modeling of the transverse 
momentum power spectrum. 

We use Equation dS]) to calculate the observable CMB 
angular power spectrum of kSZ, Ce. In the left panel of 
Figure 01 we show dCi^^Qoo/ds, which is the contribution 
to the kSZ signal at .^ = 3000 per comoving distance from 
the observer. The area under each curve gives the total 
Cf, which is shown in the right panel. The dotted line 
shows the LO calculation (Eq. ITOl) . The solid line and 
the diamonds show the Standard kSZ model (Eq. IT^ 
with the nonlinear matter densit y power spectrum com- 
puted using the HALOFIT code ([Smith et al.ll2003ll and 
the simulation (i.e., the square symbols in Figure [5]), re¬ 
spectively. We cannot calculate the kSZ contribution 
at z < 0.5 in our simulation because the contributing 
wavenumber is beyond the resolution limit of our sim¬ 
ulation (fc ~ 2 /i Mpc~^). The disagreement between 
the two at z = 0.5 is due to the disagreement between 
the HALOFIT power spectrum and our simulation at 
k > 1 h Mpc“^, which is also due to the resolution limit 
of the simulation. 

The dashed line shows the next-to-leading order con¬ 
nected moment contribution. We also estimate the con¬ 
nected term contribution in our simulation by subtract¬ 


ing the diamonds from the total transverse momentum 
power spectrum measured from the simulation. The sim¬ 
ulation at z = 1 suggests about a factor of 3 larger ef¬ 
fect than the next-to-leading order perturbation theory. 
The interpretation of the difference at z = 0.5 is difficult 
because of the disagreement between the diamond and 
the solid line; however, if we assume that the resolution 
of the simulation affects the unconnected and connected 
terms in the same way, then the simulation at z = 0.5 
suggests about a factor of three larger effect than pertur¬ 
bation theory. This is expected because the kSZ signal 
receives contributions from k > 1 h Mpc~^ where the 
next-to-leading order calculation is insufficient to cap¬ 
ture nonlinearity (see the triangles in Figured]). The 
difference decreases as z increases, as expected from de¬ 
velopment of nonlinear structure formation. Taking into 
account this extra nonlinearity above perturbation the¬ 
ory, we estimate that the connected term contribution to 
Ce at £ = 3000 is at least ten percent of the unconnected 
term. 

6 . SUMMARY AND CONCLUSION 

We have reexamined the currently popular model of 
the transverse dark matter momentum power spectrum 
(Eq. [T3|) used in the previous calcul ation of the post - 
reionization kSZ power spectrum fe.Ef.. lShaw et mMI^. 
We find that the current model reproduces the contribu¬ 
tion of the unconnected term (Eq. |9|) well. However, this 
model ignores the contribution from the connected term 
that arises in the nonlinear regime. Using both pertur¬ 
bation theory and cosmological V-body simulation, we 
show that the contribution from the connected term adds 
a significant additional power, especially at larger k at 
lower z (Figured]). This is consistent with the expecta¬ 
tion from nonlinear structure formation. 

We estimate that the contribution of the connected 
term to the kSZ angular power spectrum at £ = 3000 
is ten percent of the unconnected term. This is the 
term that needs to be added to the calculation of 
iShaw et al\ (l20I2t) . assuming that a similar correction 
is necessary for the momentum power spectrum of 
gas. In light of the current observational constraint 
(L>^£|ooo = and the post-reionization kSZ 

model of IShaw et~^ (I20I2I1 (T^^laopo^ = 2.0 
adding ten percent to the post-reionization kSZ signal 
would imply a twenty percent less k SZ signal from 
the E oR. The semi-numerical model of iBattaglia et al\ 
()20I3f) implies that such a change in the reionization 
kSZ signal would result in a twenty percent increase in 
the redshift at which the mean ionized fraction of the 
universe reached 50%, or ten percen t shorter duration of 
reionization, Az. iPark et (|2013[ ). however, show that 
reionization simulations based upon radiative transfer 
and iV-body simulation can yield a more extended dura¬ 
tion of the reionization epoch than these semi-numerical 
models predict, so this quantitative conclusion about 
the constraint from changing the kSZ upper limit from 
the EOR on the duration of reionization may need to be 
revised to take account of such extended reionization. 

Finally, we have shown that both the LO and next- 
to-leading order perturbation theory calculations give 
Pqj^ oc in the low-fc limit, independent of cosmology or 
the initial power spectrum. This behavior is consistent 



















































7 


Contributing Wavenumber (h Mpc"’) 

5.0 2.0 1.0 0.7 0,6 CMB Power Spectrum 



Fig. 4. — Angular power spectrum of the kSZ effect. (Left) dCi= 2 ,QQQ/ds^ which shows the contribution to Cf at £ = 3000 from a given 
comoving distance. The horizontal axis shows redshifts such that it is linear in the comoving distance. The top label shows the contributing 
wavenumber to the kSZ sign al, i /s{z). The discontinuity at 2 = 3 is due to the instantaneous helium reionization. The solid line shows 
the Standard kSZ model fEa. I13II with the HALOFIT density power spectrum, the das hed line shows the next-to-leading order connected 
term from perturbation theory, and t he d otted line shows the LO contribution (Eq. [Toll. We also show the simulation results at 2 = 0.5, 1, 
and 2: the diamonds show Equation {13 with the nonlinear density power spectrum from the simulation, while the crosses use the total 
simulation momentum power spectrum minus the diamonds. (Right) = £(£+ l)C£/(27r). The dash-dotted line shows the primary CMB 
temperature power spectrum, while the other lines show the same cases as in the left panel. 


with momentum conservation. 
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APPENDIX 

KSZ EFFECT FROM LONGITUDINAL MODES 
Rewriting Equation ([2l) in terms of Fourier mode of q gives 


at,.,, aTUeo /" ds 

^ 

T c J a[sY 


(Pk 


[7 ■ q(k, s)]e 


—ik- (37) 


(Al) 


Here, the momentum vector in Fourier space, q, can be decomposed into the longitudinal mode, g|| = q ■ fc, and the 
transverse component, g_L = |q — fc(q ■ k)\, to give 


at UTUefi 

T c I a(<='>2' 


/ 0^2® ^ J +cos(<(i,j-<(ip)(l-x^)^/^Qj_(k,s)]e (A2) 


where x = k ■ j. 

In this section, we shall derive the angular power spectrum of gy term. The derivation given here parallels that for 
the transverse mode given in Appendix of [Park et al\ (j2013l) . Spherical-harmonics decomposition of g|| term in the 
above equation is. 


/ AT 

d^Yr(n)—{i) 

d^vr^i) 


c 

(Pk 


ds 


d^k 


—iksx 




(A3) 
























where 


^—iksx 


fem{^) = J J ^xqii{k,s)e- 

= [ d^wr*(n) [ -J^e-'^xq\\{k,s) X 4TrJ2i-i)^jL{ks)Yl^{l)Y£^*ik). (A4) 

d J \ ) 

Here, we choose a coordinate system that k = z, thereby simplifying Yl^*(k) to Y^*{z) = Then, 

fim{kz) = V^ J -^e“^g||(k,s)^(-i)-^ji(fes) J cos6»Y^™*(7) 


IGtt^ /■ ds 




■e-"g||(k,s)^(-i)^ji(fcs) / d?Wl{^)Y?{^)Yr{l) 


(A5) 


where 6 is the azimuthal angle of 7 in the coordinate system of our choice. Also, we have used the fact that Y^{9, (j)) = 
h\f^ cos6. 

Z y TT 

The integral over 7 can be given in terms of the Clebsch-Gorden coefficient: 


I = y^^^^C'a(T,0;m,0)C«(T,0;0,0)5o,™. 


In this case, the relevant ones are 

Qi(^+1,0;0,0) = 

and, they give 


1 


2e + i 


C'n(^,0;0,0) = 0, - 1,0; 0,0) = 


2£ + V 


(A 6 ) 


(A7) 


q\\{k, s)'^(-i)^ jL{ks) x (^^(T, 0 ; 0 , 0 )C'£i(L, 0 ; 0 , 0 ) 

= H)^+V4^(2^+1) [ -^e 

J a{sy 


^ ^ I . ! ^ 


= (-f)^+^\/7r(2£ + 1) 


ds ^_^je{ks)_ 


i{sy 


ds 


h,m^o{kz)= 0 . 


(A 8 ) 


Isotropy of the universe allows us to generalize above for fimik) with any k. 
Then, we obtain the final expression for the angular power spectrum: 


C" 


/neCrT\2 1 r ds r ds' f ,, djz{ks) dj(,{ks') 

\ c ) ttJ a(s )2 J a(s')2 J ds ds' 


\JPq\\{k,z{s))Pq^^(k,z{s')). 


(A9) 


We can show that cj is much smaller than Ci in Equation Q using Pq^ and Pq^^ from the 2nd order calculation. 
For Pq ^, 2nd order is the leading order, and the expression (Eq. [TOl) and the derivation are shown in Section O For 
Pqii, the 2 nd order terms follows the 1 st order term: 


Pq\\ {k, Z) 



P^T(.k,z) + d^f 



k'\,z)p^ri^',z) 


fc/r'(fc/r' — k'fi'^ + k') 
fc'2(fe2 + k'^ _ 2fcfcV) 


(AlO) 


(|Ma a,Tid Frvll2nn2ll . Here, Pjf is the linear density power spectrum and p,' = k • k'. 

In Figure [U we show the 2nd order calculation of kSZ angular power spectrum from longitudinal and transverse 
modes. 1 st order term in the longitudinal modes dominates at ^ 100 , but diminishes rapidly in increasing (. as 

argued in iVishniad (|I987f) . At £ > 3000 where the primary CMB vanishes to the level of allowing kSZ signal to be 
measured, the longitudinal mode contribution is below the transverse mode contribution by four orders of magnitude 
or more. Based on this comparison, we do not expect the longitudinal modes to be important even to a percent level 
and we shall ignore their contribution in this work. 
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CMB Power Spectrum 



I 


Fig. 5.— Angular power spectrum of the kSZ effect from the transverse momentum field (solid) and the longitudinal momentum field 
(dot-dashed) in second order. The short-dashed and dotted lines denote contributions from the 1st and 2nd order terms in the longitudinal 
mode, respectively. The long-dashed line shows the primary CMB power spectrum. 


CONNECTED TERM IN PERTURBATION THEORY 

We take the third term in Equation ® and express it in terms of ^ = V-v, which is equivalent to say v(k) = i{)^/k)6\ 

r (fk' f dPk" . 


P — 


(27r)^ J (27r)^ 

= -{afr 


[k' • k" - (k • kO(k • k")]PM™,c(k - k', -k - k", k', k") 

P55,e.c(k-k',-k-k",k',k"). 


CO 

CO 

’k'-k"- (k-k')(k-k")' 

/ (2^)3 j 

(27r)3 

k'k" 


(Bl) 


Let us derive Pssee,c using perturbation theory. We begin with the next-to-leading order expression for the full 
expression for Pssee that includes the unconnected terms; 

(27r)^P55ee(ki,k2,k3,k4)5D(ki + k2 +k3 + k4) 

= (j(ki)5(k2)0(k3)0(k4) 

= (j^^^(ki)^P^(k 2 ) 0 ^^^(k 3 ) 0 P^(k 4 )) + cyclic (6 terms) 

^ 5 P^(ki)^P^(k 2 ) 0 ^^^(k 3 ) 0 P^(k 4 )^ + cyclic (4 terms). (B2) 

The numbers in the superscripts indicate the perturbation order. We refer to the first case as P 1221 term and the 
second case as P 3111 term. Note that above expression is not symmetric for switching one of ki and k 2 with one of 
ks and k 4 although it is symmetric for ki o k 2 and k 3 y-)- k 4 . 

One of P 1221 terms, / 5 P)(ki) 5 P)(k 2 ) 0 P^(k 3 ) 0 P)(k 4 )\, is given by 


(ki)^(P (k 2 ) 0 (^) (k3)^^^^ (k4 


" ^d^k 2 ,b / ^ / d3k3.b 


J (27r)3 J (27r)3 ^ 

C2*^(k2,aj k2,b)G2*^(k3_a, k3 b)(^D(k2 — k2^a — k2^b)fc(k3 — k3 a — k3_b) 
Jp) (ki)^(l) (k2,a)5^'^ (k2,b)5^'^ (k3.a)^'^> (ka.b)^^'^ (k4)l 
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k2,b)G2*^(k3^a, k3^b)fc(k2 — k2^a — k2,b)<^D(k3 — k3^a — ks^b) 

'(5(l)(ki)5W(k2^a)) (5('Hk2,b)^(^nk3,a)) (^(^Hk3,b)5W(k4)) 

- (^(I)(ki)5(l)(k2,a)) (5(^Hk2,b)5(l)(k3,b)) (^(^Hk3.a)5(^Hk4)) 

- (5(')(ki)^(l)(k2,a)) (5(lHk2,b)5(l)(k4)) (5(l)(k3,a)^('Hk3,b)) (vanish) 

-<Jj(i)(ki)|(i)(k 2 ,b)) (... are similar to the above 3 lines) 

- <^j(^)(ki)|(^^(k 3 ,a)^ ^^^^Hk 2 ,a)^^^Hk 2 ,b)^ (k 3 ,b)^^^Hk 4 )^ (vanish) 

- (5(l)(ki)^(l)(k3,a)) (5(^Hk2,a)^(^Hk3.b)) (k2,b)5 (^Hk4 ) ) 

■ (^(l)(ki)^(l)(k3,a)) (5(^Hk2,a)^^'Hk4)) (k2,b)^(^nk3,b) ) 

■<Jj(i)(ki)^(i)(k3,b)) (... are similar to the above 3 lines) 

(... are the unconnected terms) , 


(B3) 


where F^"\q^,c^ 2 ) = f + + f and G^"^(qi,q 2 ) = f + + | ■ The case that ^B)(ki) is 

paired with ^^^^(k 2 ,a) is equivalent to the case that it is paired with jBi(k 2 ,b), the case that ^^^^(k 2 ,a) or 0 ^^^(k 3 .a) is 
paired with ^B)(k ;2 b) or 0 ^^^(k 3 ,b), respectively, vanishes, and the case that m^i(ki) is paired with 5B)(k4) belongs 
to the unconnected moment. We proceed with the connected terms only: 



T’2*^(k2,a,k2,b)G2^^(k3 a,k3 b)l5D(k2 — k2^a ” k2^b)^D(k3 — ks a — ks b) 
'(^(ki)^(k2,a)) (5(k2,b)5(k3,a)) (k3,b)5(k4) ) 

+ (5(ki)5(k3,a)) (^(k2,a)^(k3.b)) (5(k2,b)^(k4))' 


= 4 




2,a 


(2^)3 


d^k. 


d^k. 


2,b 


3,a 


(27r)3 


d-^k; 


3,b 


T2*^(k2,a,k2,b)G2*''(k3^a,k3_b)dD(k2 — k 2 .a ~ k2.b)dD(k3 — ks^a — ks^b) 




[(27r)^P(ki)dD(ki + k2.a)(27r)^P(k2,b)fc(k2.b + k3^a)(27r)^P(k4)dD(k3,b + k4) 
+ (27r)^P(ki)(5D(ki + k3 a)(27r)^P(k2^a)^D(k3.b + k2,a)(27r)^P(k4)dD(k2,b + k4)] 


— 4(27r)^(5D(ki + k 2 + ks + k 4 ) 

[Pi"^(ki +k2,-ki)G("Hki +k2,k4)P(ki)P(ki +k2)P(k4) 

+ G(*)(ki +k3,-ki)P2(*)(ki +k3,k4)P(ki)P(ki +k3)P(k4)]. (B4) 

Similarly, 

(d(2)(ki)d(B(k2)0(2)(k3)d(l)(k4)^^ 

= 4(27r)^(5D(ki + k 2 + k3 + k4) 

[Pi"^(k2 + ki, -k2)G("^(k2 + ki,k4)P(k2)P(k2 + ki)P(k4) 

+ G^*) (k 2 + ka, -k 2 )Pi"^ (k 2 + ka, k4)P(k2)P(k2 + ka)P(k4)], 


(B5) 
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= 4(27r)^(5D(ki + k2 + ka + k4) 

[F^^\kl + k2, -ki)G^"^(ki + k2,k3)P(ki)P(ki + k2)P(k3) 

+ G^")(ki +k4,-ki)P2(")(ki +k4,k3)P(ki)P(ki +k4)P(k3)], 
id 

(p)(ki)5(i)(k2)0(i)(k3)0('Hk4))^ 

= 4(27r)^(5j3 (ki + ]^2 + ks + k 4 ) 

[Pi"Hk 2 +ki,-k 2 )G("Hk 2 +ki,k3)P(k2)P(k2 +ki)P(k3) 

+ G^") (k 2 + k4, -k 2 )Pi"^ (k 2 + k4, k3)P(k2)P(k2 + k4)P(k3)]. 

The following two terms can be expanded in the same way as well, but have slightly different forms: 

(P)(ki)5(2)(k2)0(l)(k3)0(l)(k4))^ 

= 4(27r)^(5D (ki + k 2 + ks + k 4 ) 

[Pi"^(k3 +ki,-k3)Pi"^(k3 +ki,k4)P(k3)P(k3 +ki)P(k4) 

+ Pi")(k3 + k 2 , -k3)Pi"^(k3 + k 2 , k4)P(k3)P(k2 + k3)P(k4)], 


( (ki)5(l) (k2)d(2) (k3)0(") (k4))^ 

= 4(27r)^5D (ki + k2 H- ks + k4) 

[G^"^(ki +k3,-ki)G("^(ki +k3,k2)P(ki)P(ki +k3)P(k2) 

+ G^*)(ki + k4, -ki)G^*)(ki + k4, k2)P(ki)P(ki + k4)P(k2)]. 

For P 3111 terms, we have 

(p)(ki)5(i)(k2)0'i)(k3)0(')(k4)) 

r dk^„ f dk^, f , 

= j j J dfc3cP3(ki.a,ki,b,ki.c)<5r3(ki-ki,a-ki,b-ki,c) 

( (ki,a)5(i) (knb)^('> (ki.c)^(^’ (k2)5(i) (k3)5(l) (k4)) 

r dk^ „ f dkF r ^ 

= J J ^ J P3(ki,a,ki,b,ki,c)<5i3(ki-ki,a-ki,b-ki,c) 

(ki,a)5(i) (knb)^^^^ (ki,c)^(^^ (k2)5(i) (k3)5(l) (k4)) 

= 6(27r)^P3 ^(—k2, —ka, —k4)P(k2)P(k3)P(k4)(5£i(ki + k2 + ka + k4). 


(B 6 ) 


(B7) 


(B 8 ) 


(B9) 


(BIO) 


(j(l)(ki)5(3)(k2)0(l)(k3)0(l)(k4)) 

= 6(27r)^P3 ^(—ki, —ka, —k4)P(ki)P(k3)P(k4)(5£)(ki + k2 + ka + k4), (Bll) 


(j(l)(ki)5(l)(k2)0(3)(k3)0(l)(k4)) 

= 6(27r)3G^"^(-ki, -k2, -k4)P(ki)P(k2)P(k4)<5D(ki + k2 + ka + k4), (B12) 

(j(l)(ki)5(l)(k2)0(l)(k3)0(3)(k4)) 

= 6(27r)3G("^(-ki, -k 2 , -k3)P(ki)P(k2)P(k3)5D(ki + k 2 + ka + k4). 


(B13) 
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where the kernels, F 3 and G 3 , are given by the following recursion relations. 


F„(qi,...,q„)= V , [( 2 n+ l)a(qi + ... +qm,qm+l + ••• + qn)^’n-m(qm+l, ■■■, Qn) 

(2n + 3)(n- 1) 

-t-2/3(qi “t" ... -j- qm; qm+l 'b ... ”f~ qn)G7i_77j(qm-|-i,..., qn)] : 


G„(qi 


N Gm(qi,..., qm) ro ^ , , , , 

qn) — / , fn I o\/ 1 \ [^'^(qi + ■■• + qmj qm+l + ... + qn)^n—m 

(2n + 3)(n- 1) 

+2n/3(qi + ... + q^, qm+i + + qn)G n—m (qm+l,---,qn)] , 


(qm+l; qn) 


(B15) 


(s') (s') 

and Fg ' and Gg ' denote the symmetrized kernels of Fg and G 3 . _ 

Combining results above and substituting them into Equation (IBI 6 II . we obtain the expression for the transverse 
momentum power spectrum from the connected terms: 


pNLO 


(Pk' 


r Pk" 

’k'-k"- (k-kO(k-k")' 

1 (27r)3 

fc'/c" 


Fg(^)(k + k",-k', 


-k")F(k + k")F(k')F(k") 


+ Fg^'*^ (k - k', k', k")F(k - k')F(k')F(k") 

+ G^'*^ (k - k', -k - k", k")F(k - k')F(k + k")F(k") 
+ G^'*^ (k - k', -k - k", k')F(k - k')F(k + k")F(k')] 


+ 4[F2^'’^(k - k', k' + k")G^"^(k' + k", -k")F(k - k')F(k' + k")F(k") 

+ F2^'*^ (k - k', k' + k")G^"^ (-k' - k", k')F(k - k')F(k' + k")F(k') 

+ G^"^ (k, -k + k')F2^"’ (k, k")F(k - k')F(k)F(k") 

+ G^"^ (k, -k + k')G^"^(k, -k - k")F(k - k')F(k)F(k + k") 

+ (k - k' + k", -k + k')F2^'’^ (k - k' + k", k')F(k - k')F(k - k' + k")F(k') 

+ G^'’^ (k - k' + k", -k + kOG^'^^ (k - k' + k", -k - k")F(k - k')F(k - k' + k")F(k + k") 
+ F2^'’^(-k' - k",k + k")G^"^(-k' - k", k")F(k + k")F(k' + k")F(k") 

+ F2^'’^(-k' - k",k + k")G^"Vk' - k", k')F(k + k")F(k' + k")F(k') 

+ G^'’^ (-k - k" + k', k + k")F2^"^(-k - k" + k', k")F(k + k")F(k - k' + k")F(k") 

+ G^"^ (-k, k + k")F2^'’) (-k, k')F(k + k")F(k)F(k') 

+ F!f'> (k, -\^')F!f'> (k, k")F(k)F(k')F(k") 


+ F2^'*^(k - k' + k", k')F2^"^(k - k' + k", 


-k")F(k')F(k - k' + k")F(k") 


} 


(B16) 


UNCONNECTED TERM IN PERTURBATION THEORY 

( 2 ) (2") (2") 

For the unconnected term, we begin by substituting the next-to-leading order power spectra of Pgg' and Pgg' 
in Equation (121). Then, perturbation theory gives 


Ps?(M)= j |^Fi'-^(k-k',k')2F(|k-k'|)F(F) 

+ 2 J -0^Ft\k,-k',k')P{k)P{k'), 

Ps!\^)= J -0^Ft\k-k',k')Gi^\k-k\k')Pi\k-k'\)Pik') 
+ J ^^i'^(k,-k',k')F(fc)F(fc') 

















13 


Pel\^)= I -0^Gi^\k-k',k'fP{\k-k'\)Pik') 

+ 2 J (k, -k', k')P{k)P{k'). (Cl) 

Substituting Pss = PggK Pse = Ps^ and Pgg = Pg^ to Equation ([5]) and dropping higher order terms like p( 2 )p( 2 )^ 
we obtain the expression for the transverse momentum power spectrum from the unconnected terms: 


Pq!^^c{k, z) = j ^^(1 - 


k'^ 


P(|k - k'\)P^f{k^) + P^'t>i\k - k'\)Pik') 


( 2 )/ 


|k-k'|2 L 


P(|k - k'DP^fik') + P^'^g>i\k - k'\)P{k') 


(2), 


(C2) 
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